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Abstract 

Within the framework of the improved isospin dependent quantum molecular dynamics 

(ImlQMD) model, the pion emission in heavy-ion collisions in the region 1 A GeV is investi- 

\Q ' gated systematically, in which the pion is considered to be mainly produced by the decay of 
OS '■ 

. resonances A(1232) and N*(1440). The in-medium dependence and Coulomb effects of the 

I pion production are included in the calculation. Total pion multiplicity and vr /vr^ yields are 

■ calculated for the reaction ^^'^Au-l-^^'^Au in central collisions for selected Skyrme parameters 



0\ 

o 



SkP, SLy6, Ska, SIII and compared them with the measured data by the FOPI collaboration. 



^ : PACS: 25.75.-q, 13.75. Gx, 25.80.Ls 
P.. 

Heavy-ion collisions at intermediate energies play a significant role to extract the information 
of the nuclear equation of state (EoS) under extreme conditions, i.e., at high densities and high 
temperature. Besides nucleonic observables such as rapidity distribution and flow, also mesons 
emitted from the reaction zone can be probes of the hot and dense nuclear matter, that are also the 
interest physics at the Cooling Storage Ring (CSR) energies in Lanzhou. ^ The emission of pion in 
heavy-ion collisions in the region 1 A GeV is especially sensitive as probes of isospin asymmetric 
EoS at supra-saturation densities.^ Spectra of the pion emission in heavy-ion collisions have 
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been measured by the Kaos and FOPI collaborations and analyzed systematically by the present 
theoretical transport models.^ ^ A comparison of the various transport approaches was made 
in Ref. [5J. The present theoretical models overpredict the total pion multiplicity if using free 
nucleon-nucleon (NN) cross sections below 2 A GeV region compared with experimental data. 
Further investigations of the pion emissions in the 1 A GeV region are still necessary by improving 
transport models or developing some new approaches. The improved isospin- dependent quantum 
molecular dynamics model has been successfully applied to treat fusion dynamics and reaction 
mechanism of two colliding nuclei near Coulomb barrier. 1^ ^ 1^ To investigate the pion emission, 
we further include the inelastic channels in nucleon-nucleon collisions in the ImlQMD model. 

In the ImlQMD model, the time evolutions of the baryons and pions in the system under the 
self-consistently generated mean-field are governed by Hamilton's equations of motion, which read 
as 

OH .OH 

Here we omit the shell correction part in the Hamiltonian H as described in Ref. [7|. The 
Hamiltonian of baryons consists of the relativistic energy, the effective interaction potential and 
the momentum dependent part as follows: 



HB=Y.yPi+ ^^ + Uint + U^om ■ (2) 
i 

Here the Pj and rrii represent the momentum and the mass of the baryons. 

The effective interaction potential is composed of the Coulomb interaction and the local inter- 
action 

Uint = Ucoul + Uloc- (3) 

The Coulomb interaction potential is written as 

Ucoui = ^ E ^erf{r,,/VAL) (4) 

where the ej is the charged number including protons and charged resonances. The Vij = \ri — rj\ 
is the relative distance of two charged particles. The local interaction potential is derived directly 
from the Skyrme energy- density functional and expressed as 

Uloc = J VioMr))dr. (5) 
The local potential energy- density functional reads [7] 

T/ f ^\ ^ , P^^^ I 9sur n2 I fl'surrvT/ m2 , gym 2 c-2 , 8/3/ 5/3 /^n 

yiocip) = IT — + 1 , ,7 + ^(^^) + -^y^^P^ - PP)\ + -^P ^ +9rp' /Po , (6) 

^ Po -L + 7 Po '^Po '^Po ^Po 
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where the p is the baryon density and the 5 = (pn—Pp) / [Pn + Pp) is the isospin asymmetry with the 
proton density pp and the neutron density pn ■ The momentum dependent part in the Hamiltonian 
is expressed as 

Umom=^- E ^Ne(p,-p,.)' + l)F, (7) 



with 

r f.._,^ 



AL 

Here the L denotes the square of the pocket-wave width, which is dependent on the size of the 
nucleus. 

In Table 1 we list the ImlQMD parameters related to several typical Skyrme forces after 
including the momentum dependent interaction. The parameters a, /3, 7, Qr, gsur, Qlun ^ 
e are related to the Skyrme parameters ^0,^1,^2,^3 and xq, xi, X2, X3, and determined in order to 
reproduce the binding energy (£'^=-16 MeV) of symmetric nuclear matter at saturation density 
for a given incompressibility as well as the correct momentum dependence of the real part of the 
proton-nucleus optical potential. In the following calculation we take the Skyrme parameter SLy6, 
which can give the good properties from finite nucleus to neutron star [9J. 

Analogously to baryons, the Hamiltonian of pions is represented as 



i=l 



H. = Y.( VP? + < + ' (9) 



where the and represent the momentum and the mass of the pions. The Coulomb interaction 
is given by 



yCoul ^ ^ (^Q) 

where the N^^ and Nb is the total number of pions and baryons including charged resonances. 
Thus, the pion propagation in the whole stage is guided essentially by the Coulomb effect. The 
in-medium pion potential in the mean field is not considered in the model. However, the inclusion 
of the pion optical potential based on the perturbation expansion of the A-hole model gives 
negligible influence on the transverse momentum distribution. 

The pion is created by the decay of the resonances A (1232) and N*(1440) which are produced 
in inelastic NN scattering. The direct pion production cross section is very small in the considered 
energies and not included in the model. The reaction channels are given as follows: 

NN^N/\, NN^NN\ NN ^ AA, 

A^Nn, N*^N7i. (11) 



The cross section of each channel to produce resonances are taken the values calculated with the 
one-boson exchange model. ^ Transport models overpredicted the total pion production with the 
free cross section. In the ImlQMD model, we use the free elastic cross section and the in-medium 
inelastic cross section which is given by cr^Zfdiwm — iji^y'^Jree^^^^ with the free baryon-baryon 
(BB) inelastic cross section aj^f^^^^^ and the reduced effective mass /i^^ (free mass Hbb)- The 
experimental data of total elastic and inelastic cross sections^I^ are parameterized in the ImlQMD 
model as shown in Fig.l. 

In the 1 A GeV region, there are mostly A resonances which disintegrate into a tt and a 
nucleon, however, the A^* yet gives considerable contribution to the high energetic pion yield. The 
energy and momentum dependent decay width is used in the ImlQMD model and expressed as 

(1 + a2|p|^)(a3 + IpI^) 

which originates from the p-wave resonances. The p is the momentum of the created pion (in 
GeV/c) in the resonance rest frame. The values ai =22.48 (17.22), 02 =39.69 and 03 =0.04 (0.09) 
are used for the A (A^*) with bare decay width Tq =0.12 GeV (0.2 GeV). [13] in Fig.2 we show a 
comparison of the time evolution of the vr, A and A^* production in the reaction ^^''Au+^^^Au for 
head on collisions at 1 A GeV for two cases of the bare decay and energy dependent decay widths. 
Both methods almost give the same yield of the pion production. In the following calculation, 
we use the energy and momentum dependent decay width. In Fig.3 we give the multiplicity of 
produced pion as a function of the impact parameter for the same system at 1 A GeV energy. The 
numbers of produced vr^, vr'^and 7r+ are reduced with increasing the impact parameter because of 
the decrease of the participants of the 'fire ball' formed in the heavy-ion collisions. 

The emission of the produced pion is sensitive to the incident energy owing to the size of the 
compressed nuclear matter. We calculated the transverse momentum distribution of vr", vr^and vr^ 
in central ^^''Au+^^^Au collisions at different incident energies as shown in Fig. 4. The larger and 
wider distributions were found at the higher incident energies due to the larger participant numbers 
of the collision nucleons. The high energy pions originate from the early phase and the decay of the 
N*(1440) resonance also plays a significant role.^^^ In Fig. 5 we compare the total pion number and 
the 7r~/7r+ ratio with the FOPI data in central ^^^Au-|-^^^Au coUisions^Sl for the Skyrme parameters 
SkP, SLy6, Ska and SIII which correspond to the different incompressibility modulus as listed in 
table 1. The calculated value of the total pion number is related to the incompressibility modulus 
Kao and the effective mass in nuclear medium. Over the whole domain, the force SLy6 is nice 
and can reproduce the experimental data. But the parameter slight overpredicts the total pion 



multiplicity at lower incident energies and underestimates the value at higher incident energies 
if using the above in-medium inelastic cross section. The in-medium elastic and inelastic cross 
sections are still open problems in transport model calculations, which should be calculated by 
microscopic many-body models and then parameterized to add into transport models. The 7r~/7r+ 
ratio is interest for extracting the high density behavior of the symmetry energy per nucleon.^^ 
Using the isobar model, one gets the ratio 7r^/7r+=1.95 for pions from the A resonance, and 
7r^/7r+=1.7 from the A^* for the system ^^^Au+^^''Au.EIS These relations are globally valid, i.e. 
independent of the pion energy. The observed energy dependence of the 7r^/7r+ ratio is due 
to the influence of the Coulomb force and the symmetry energy interaction. The tt^/tt^ ratio 
is sensitive to the stiffness of the symmetry energy at the lower incident energies. Recently, a 
soft nuclear symmetry energy at supra- saturation densities was pointed out by fitting the FOPI 
data with the IBUU04 model. EISl in the ImlQMD model, we only consider the linear dependence 
of the symmetry energy term on the baryon density as shown in Eq.(6). The inclusion of the 
density-dependent symmetry energy in the ImlQMD model is in progress. 

In summary, the pion production in heavy-ion collisions in the region 1 A GeV for the reaction 
^^^Au-|-^^^ is investigated systematically by using the ImlQMD model. The distribution of the 
transverse momentum is calculated at different incident energies. The total number of produced 
pion and the 7r^/7r+ ratio are calculated in central collisions for selected Skyrme parameters SkP, 
SLy6, Ska, SIII and compared them with the FOPI data. Deviations from the simple isobar model 
originate from the Coulomb and the symmetry interactions. The Tr'/Tr"*" ratio is sensitive to the 
stiffness of the sjnnmetry energy at the lower incident energies that may be further investigated 
at the CSR energies. 

We would like to thank Prof. Lie- Wen Chen, Prof. Wei Zuo and Dr. Gao-Chan Yong for 
fruitful discussions. 
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Table 1: ImlQMD parameters and properties of symmetric nuclear matter for Skyrme effective 
interactions after the inclusion of the momentum dependent interaction with parameters 5=1.57 
MeV and e=500 cVGeV^ 



Parameters 


SkM* 


Ska 


SIII 


SVI 


SkP 


RATP 


SLy6 


a (MeV) 


-325.1 


-179.3 


-128.1 


-123.0 


-357.7 


-250.3 


-296.7 


(3 (MeV) 


238.3 


71.9 


42.2 


51.6 


286.3 


149.6 


199.3 


7 


1.14 


1.35 


2.14 


2.14 


1.15 


1.19 


1.14 


w(MeV fm2) 


21.8 


26.5 


18.3 


14.1 


19.5 


25.6 


22.9 


^r(MeV fm^) 


-5.5 


-7.9 


-4.9 


-3.0 


-11.3 


0.0 


-2.7 


Qr (MeV) 


5.9 


13.9 


6.4 


1.1 


0.0 


11.0 


9.9 


Csym (MeV) 


30.1 


33.0 


28.2 


27.0 


30.9 


29.3 


32.0 


Poo (fm-3) 


0.16 


0.155 


0.145 


0.144 


0.162 


0.16 


0.16 


m*^/m 


0.639 


0.51 


0.62 


0.73 


0.77 


0.56 


0.57 


(MeV) 


215 


262 


353 


366 


200 


239 


230 



proton-proton/neutron-neutron 
• total (exp.) 
o elastic (exp.) 
■ inelastic (exp.' 

total (ImlQMD) 

-elastic (ImlQMD) 

inelastic (ImlQMD) 




10'' 10° 10 

E..(GeV) 



E 



proton-neutron 
• total (exp.) 
o elastic (exp.) 
■ inelastic (exp.) 

total (ImlQMD) 

elastic (ImlQMD) 

inelastic (ImlQMD) 




Figure 1: Comparison of nucleon-nucleon cross sections parameterized in ImlQMD and the ex- 
perimental data.f 
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Figure 2: Production of pion, delta and N* for head-on collisions in the reaction ^^^Au+^^^Au at 
1 A GeV as functions of evolution time with the energy dependent decay width (left panel) and 
fixed width (right panel). 




Figure 3: Final multiplicities of tt , 7r° and tt"*" as a function of impact parameter for head-on 
collisions in the reaction ^^^Au-|-^^''Au at 1 A GeV. 
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Figure 4: Final transverse momentum distribution in central ^^^Au+^^^Au collisions at different 
incident energies. 
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Figure 5: Calculated excitation functions of the total pion multiplicity (left panel) and the ratio 
7r~/7r+ (right panel) in central ^^^Au+^^^Au collisions and compared with the FOPI data.l^l 



